Beyond scaling and locality in turbulence 
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An analytic perturbation theory is suggested in order to find finite-size corrections to the scaling 
power laws. In the frame of this theory it is shown that the first order finite-size correction to 
the scaling power laws has following form S(r) = cr a °[\n(r /rf)] ai , where 77 is a finite-size scale (in 
particular for turbulence, it can be the Kolmogorov dissipation scale). Using data of laboratory 
experiments and numerical simulations it is shown shown that a degenerate case with cto = can 
describe turbulence statistics in the near-dissipation range r > rj, where the ordinary (power-law) 
scaling does not apply. For moderate Reynolds numbers the degenerate scaling range covers al- 
■ most the entire range of scales of velocity structure functions (the log-corrections apply to finite 

Reynolds number). Interplay between local and non-local regimes has been considered as a possi- 
, ble hydrodynamic mechanism providing the basis for the degenerate scaling of structure functions 

^v^j . and extended self-similarity. These results have been also expanded on passive scalar mixing in 

turbulence. Overlapping phenomenon between local and non-local regimes and a relation between 
[ position of maximum of the generalized energy input rate and the actual crossover scale between 

these regimes are briefly discussed. 

PACS numbers: 47.27.-i, 47.27.Gs 
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INTRODUCTION 




Scaling and related power laws are widely used in physics. In real situations, however, scaling holds only approx- 
imately. As a consequence, the corresponding scaling power laws hold approximately as well. On the other Jhand, 
discovery of the extended self-similarity in turbulence [H,[3] and in the critical phenomena (see for a review 
shows that universal laws similar to the scaling ones can go far beyond the scaling itself, both for velocity [1| 
and for different fields convected by turbulence [5], 6]. The question is: does the extended self-similarity (ESS) can 
be theoretically explored in the frames of general scaling ideas or one needs in a completely new frames to explain 
the ESS? In turbulence this problem is related to the problem of the so-called near-dissipation range. Theoretical 
£ — approach to description of the scaling (the so-called inertial) range of scales in turbulence requires that the space scales 
r 3> rj, where 77 is (Kolmogorov) dissipation or molecular viscous scale 0,@|. On the other hand, the near-dissipation 
range of scales, for which r > rj has a more complex, presumably non-scaling dynamics (see, for instance, Refs. |3|-[T5T| 
and the references cited there) Indeed, for moderate Reynolds numbers, the near-dissipation range can span most of 
the available range of scales Different perturbation theories, which start from scaling as a leading mode, can be 
developed. Effectiveness of such theories is usually checked by comparison with experiments. In the present paper 
we will explore a generic analytic expansion of scaling in direction of the dissipation scale rj. Effectiveness of the 
suggested perturbation theory has been shown for the turbulent flows with moderate Reynolds numbers. Interplay 
between local and non-local regimes has been considered as a possible hydrodynamic mechanism providing the basis 
, for the degenerate scaling of structure functions and extended self-similarity in turbulent flows. 
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II. ANALYTIC PERTURBATIONS TO SCALING 

Let us consider a dimensional function S(r) of a dimensional argument r. And let us construct a dimensionless 
function of the same argument 

«(r) = ^§ • (1) 
r dr 

If for L > r > i( we have no relevant fixed scale (scaling situation), then for these values of r the function a(r) must 
be independent on r, i.e. a(r) ~ const for L > r > ?). Solution of equation (1) with constant a can be readily found 
as 

S(r) * cr a (2) 



where c is a dimension constant. This is the well-known power law corresponding to the scaling situations. 
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FIG. 1: The second order structure function ^(r/r;) against r/rj. The experimental data (R\ = 206, 487) [171 ] are shown as 
circles. The solid curves are the best fit of (7) to the data, corresponding to the degenerated scaling. 



Let us now consider an analytic theory, which allows us to find generic corrections of all orders to the approximate 
power law, related to the fixed small scale rj. In the non-scaling situation let us denote 



/ = HS/A), 



x = \yl(t jr\) 



where A and r\ are dimensional constants used for normalization. 
In these variables, equation (1) can be rewritten as 

df , , 
— = a{x) 
ax 



(3) 



(4). 



In the non-scaling situation a; is a dimensionless variable, hence the dimensionless function a(x) can be non-constant. 
Since the 'pure' scaling corresponds to i > 1 we will use an analytic expansion in power series 



a(x) = «o H 

x 



(5) 



where a n are dimensionless constants. After substitution of the analytic expansion (5) into Eq. (4) the zeroth order 
approximation gives the power law (2) with a = a . First order analytic approximation, when one takes only the two 
first terms in the analytic expansion (5), gives 



S{r) = cr ao [\n(r/r])} ai . 



(6) 



This is a generic analytic first order approximation to the scaling power laws provided by the perturbation theory 
suggested above. 

Corrections of the higher orders can be readily found in this perturbation theory. 

For degenerate case with ao = the second term in the expansion (5) becomes the leading term 



S(r) ~ [\n(r/r])} ai 



(7) 



It was already mentioned that there are many different ways to develop perturbation theory to scaling. The idea 
that just hi(r/r)) is an appropriate parameter for the finite-size corrections to scaling in turbulence was suggested in 
[HI]. On a physical level this was related in [l5| to the instabilities of the thin vortex tubes (or filaments), which are 
the most prominent hydrodynamical elements of turbulent flows. These instabilities usually appear as kinks readily 
transforming into wave packets propagating along the filaments. One can consider such wave packet with scale r 
and make use of the "localized induction" approximation to the Biot-Savart formula [l6l ]. In this approximation, 
contributions of portions at distances greater than ~ r S> rj to the velocity fluctuation in the immediate vicinity of 
any given point on the filament are neglected. The vortex core radius, which usually associated with rj, provides a 
cutoff from below. The vortex filament dynamics in this approximation is described by the equation 



dt 4tt { ln (77) } 7t> 
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FIG. 2: The second order structure function S2{r/rj) against r/r). The DNS data pH ] (R\ = 460) are shown as circles. The 
solid curve is the best fit of (7) corresponding to the degenerated scaling. 

where X is the position vector of a point on the filament, T is the vortex strength, 7 is the local curvature, and b is 
the unit binormal vector of the filament. One can see that in this approximation the dependence on r is completely 
determined by the logarithmic term ln(r/r/) in the dynamic equation. This pure hydrodynamic consideration was 
suggested in }15j | as a basis for using just ln(r/77) in finite-size corrections to scaling in turbulence (see also Section IV 
below) . 

III. TURBULENCE 

The longitudinal structure function of order p for the velocity field 0] 

S p (r) = (\Au\n (8) 

(Am = (u(r + x) — u(r)) • r/r) calculated for p = 2, using the data obtained in a wind tunnel at R\ = 206 and 487 
[17j . is shown in Fig. 1. Here R\ is the so-called Taylor microscale Reynolds number. The experiment with a flow, 
which was a combination of the wake and homogeneous turbulence behind a grid, is described in Ref. (l7j . We invoke 
Taylor's hypothesis 0] to equate temporal statistics to spatial statistics. 

The solid curves in these figures are the best fit by equation (7) corresponding to the degenerate scaling. Figure 2 
shows S2 (r) calculated using data from a high- resolution direct numerical simulation of homogeneous steady three- 
dimensional turbulence (la |. corresponding to 1024 3 grid points and R\ — 460. The solid curve in this figure also 
corresponds to the best fit by equation (7). As far as we know Eq. (7) was suggested as an empirical approximation of 
the structure functions for the first time in our earlier paper [l2l ] . Application of these results to correlation functions 
Q is in good agreement with the structure functions analysis (cf. for instance, Figs. 1 and 3). 

Figure 4 shows the structure functions of different orders p for the wind-tunnel data (R\ = 487) and the degenerate 
scaling (7) is shown in the figure as the solid curves. One can see that the degenerate scaling requires two fitting 
constants just as the ordinary scaling does, and, in the examples discussed above, the range of scales covered by the 
degenerate scaling is about two decades. 

IV. A POSSIBLE HYDRODYNAMIC SCENARIO 

In order to understand a possible hydrodynamic mechanism of the phenomenon described above let us recall that 
in isotropic turbulence a complete separation of local and non-local interactions is possible in principle. It was shown 
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FIG. 3: Correlation function for velocity fluctuations against logr/rj for Rx = 487 (circles) 17j. The solid curve is drawn in 
the figure to indicate agreement with the degenerate scaling shown in Fig. 1. 

by Kadomtsev [20] that this separation plays a crucial role for the local Kolmogorov's cascade regime with scaling 
energy spectrum 

E(k) ~ e 2 / 3 fc- 5 / 3 (9) 

where e is the average of the energy dissipation rate, e, k ~ 1/r is the wave-number. This separation should be 
effective for the both ends. That is, if there exists a solution with the local scaling (9) as an asymptote, then there 
should also exist a solution with the non-local scaling asymptote. Of course, the two solutions with these asymptotes 
should be alternatively stable (unstable) in different regions of scales. It is expected, that the local (Kolmogorov's) 
solution is stable (i.e. statistically dominating) in inertial range (that means instability of the non-local solution in 
this range of scales). 

Roughly speaking, in non-local solution for small scales r only non-local interactions with large scales I (1 > r/i) 
are dynamically significant and the non-local interactions is determined by large scale strain/shear. This means that 
one should add to the energy flux e-parameter (which is a governing parameter for the both solutions) an additional 
parameter such as the strain s for the non-local solution. As far as we know it was noted for the first time by 
Nazarenko and Laval (2l| that dimensional considerations applied to the non-local asymptotic result in the power-law 
energy spectrum 

E(k) ~ - k' 1 (10) 
s 

both for two- and three-dimensional cases. Linear dependence of the spectrum (10) on s is determined by the linear 
nature of equations corresponding to the non-local asymptotic that together with the dimensional considerations 
results in (10) [2l[. Interesting numerical simulations were performed in [22j ■ In these simulations local and non-local 
interactions have been alternatively removed. For the first case a tendency toward a spectrum hatter than '-5/3' is 
observed near and beyond the separating scale (beyond which local interactions are ignored), that supports Eq. (10) 
(see also below). 

Following to the perturbation theory suggested in Section II both local and non-local regimes can be corrected. 
The first order correction is 

E(k) ~fc- 5 / 3 [ln(fc d /fc)] 7 (11) 

and 



E(k) ~ k-^Hka/k)}? 



(12) 
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FIG. 4: As in Fig. 1 (R\ — 487) but for different orders of the structure function. 



(where k ( i = 1 /rf) for the local and non-local regimes respectively (7 and fj are dimensionless constants) . In the case 
when /J > 7 the fist order finite-size correction in the non-local regime becomes substantial much 'earlier' (i.e. for 
smaller k) than for the local (Kolmorov's) regime (see Eq. (5)). This can result in viscous stabilization of the non- 
local regime (cf [Hj]) and, as a consequence, in the so-called 'exchange of stability' phenomenon at certain k c = l/r c . 
That is, for k < k c the Kolmogorov's regime is stable and the non-local regime is unstable, whereas for k > k c the 
Kolmogorov's regime is unstable and the non-local regime is stable. For this scenario, at k = k c the Kolmogorov's 
regime is still asymptotically scale- invariant (i.e. Eq. (9) gives an adequate approximation for this regime), while for 
the non-local regime the first order correction is substantial (i.e. Eq. (12) should be used at k = k c for the non-local 
regime). In this scenario the Kolmogorov's regime plays significant role in the viscous stabilization of the non-local 
regime for k > k c (cf [22]), but for these k the non-local regime becomes statistically dominating instead of the 
Kolmogorov's one. Figure 5 shows a spectrum measured in nearly isotropic turbulence downstream of an active grid 
at R\ = 716 (the data are reported in [lj|). The solid curve is drawn in the figure to indicate correspondence of the 
data to the equation (12) (non-local regime), while the dashed straight line indicates the Kolmogorov's "-5/3" law to 
show possible interplay between the Kolmogorov's and non-local regimes at this Reynolds number. 
If one tries to estimate the second order structure function S^r) as [7| 

S 2 (r) ~ / h E(k)dk (13) 

Jl/r 

for r < r c using Eq. (12), one obtains 

S 2 {r) ~ [m(r/r?)]^ (14) 

with £2 = P + 1 (cf Section 3). For instance, for the data shown in Fig. 5 we obtain £2 ~ 3.9. Comparing figures 1,2 
and 5 one can see that a is monotonically increased with Re\. 
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FIG. 5: Longitudinal energy spectrum against krj in log-log scales. The data (circles) correspond to a nearly isotropic (active-) 
grid turbulence R\ = 716 [XSI | . The solid curve is drawn in the figure to indicate agreement with the Eq. (12), while the dashed 
straight line indicates the "-5/3" Kolmogorov's law. 



Let us generalize (13) introducing an effective spectrum E n (k) 

ri/v 

S„(r) ~ / E n (k)dk (15) 

Jl/r 

(i?2(fc) = E(k)). Then using the dimensional considerations we obtain for the non-local regime 

E n (k) ~ (ij k- 1 (16) 

The first order correction is 

/ c \ n /2 

(~) k-'Mkd/k)} " (17) 

Substituting (17) into (15) one obtains 
with ( n — p n + 1. 

Thus the non-local regime can provide hydrodynamic basis for the degenerate scaling of the structure functions. 



V. EXTENDED SELF-SIMILARITY 

In our previous paper [l2T | we discussed possible relation between ESS and structure functions given by Eq. (7). 
Now we can elaborate this relation with more details. While in the near-dissipation range we observe the degenerate 
scaling an ordinary scaling (presumably Kolmogorov's o ne) has been certainly observed in the inertial range of scales 
for very large Reynolds numbers (see, for instance, 0]j[13|)- Therefore, one can expect the existence of a crossover 
scale r c from the ordinary scaling to the degenerate one (see previous section). At the scale r — r c we can use a 
continuity condition between ordinary scaling (2) and degenerate scaling (7) 



erf = A[\n{r c /rj]^ 



(19) 
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FIG. 6: Normalized exponents Cp/C3 against p for R\ = 487 (circles). The exponents were obtained using the data shown in 
Figs. 1,4. Crosses are the ESS exponents obtained for the atmospheric turbulence data at R\ = 10340 [231 ]. 



for the structure functions of all orders (c and A are some constants) . Let us denote in general case 

S p (r) = A' p A; S p (r) = A p [\n{r / 



(20) 



in inertial and near-dissipation ranges respectively. 

Then condition of compatibility of the continuity equations (19) for different order (m and n) structure functions 



A' m r£> = A m [ln(r c /r,)}^; A' n r& = A n [ln(r c / V )]^ 



(21) 



can be of two kinds. 

The first kind is applicable for arbitrary value of r c and has the form 



c 



A> \Cn 

A, 



A, \ C« 

Ay. 



(22) 



while another kind is applicable for a fixed value of r c only. Below we will be interested in the first kind of the 
compatibility condition. This condition, in particular, gives a power-law relation between moments of different orders 



with the same exponent 



S m (r) ~ S£»"(r), 



a _ Cm _ Cm 
Pm.n — t. 



(23) 



(24) 



for both inertial and near-dissipation ranges (see Eq. (22)). The last phenomenon was previously observed in different 
turbulent flows and was called Extended Self-Similarity (ESS) [l|,[3] .[EH]. It is important to note that at ESS not only 
the scalings (23) are the same in both regimes inertial and dissipative, but also the prefactors are the same. Figure 
6 shows as circles the exponents £ p obtained from Figs. 1,4 and normalized by the exponent £3- For comparison we 
show by crosses in this figure the normalized exponents obtained using ESS in the atmospheric turbulence for large 
i? A = 10340 [H|. 

It should be noted that a high-resolution numerical simulation reported in a recent paper [131 ] for low- Reynolds- 
number (R\ « 10 — 60) flows shows no hint of scaling-like behavior of the velocity increments even when ESS is applied. 
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That should return us to the Section II and recall us that the degenerate scaling in the form of Eq. (7) is only a 
first non-trivial term in the perturbation theory. Moreover, the question of convergence of the perturbation theory 
used in Section II become significant at the low-Reynolds-number, when the dissipation scale is effectively too close to 
the scales under consideration. Though, the authors of Ref. [l3[ claim that: 'the DNS scaling exponents of velocity 
gradients agree well with those deduced, using a recent theory of anomalous scaling, from the scaling exponents of 
the longitudinal structure functions at infinitely high Reynolds numbers. This suggests that the asymptotic state of 
turbulence is attained for the velocity gradients at far lower Reynolds numbers than those required for the inertial 
range to appear'. This circle of problems can be also related to the possibility of fluctuations of the dissipation scale 
[HI ! H3 ■ For the cases when these fluctuations are significant a set of the expansions might be needed or (alternatively, 
for moderate Reynolds numbers) calculation should involve averaging over the fluctuating cut-off rj. In the last case the 
experimental and numerical data suggest that the average value of rj should not deviate strongly from the Kolmogorov 
dissipation scale (at least for the moderate values of Reynolds number). 

Though, the problem with the small Reynolds numbers can be deeper. We do not know whether the non-local 
regime appears at the same Reynolds number as the local one. The energy balance of the non-local regime is more 
complex than that of the local one. Actually the non-local regime needs in the local one (even statistically unstable) 
for this balance (cf 22]). Therefore, critical Reynolds number for appearance of the (stable) non-local regime can be 
larger than that for the local regime. 



VI. GENERALIZED ENERGY INPUT RATE 



Figures 2 and 5 show that for sufficiently large Reynolds numbers, providing a visible inertial interval, there is an 
overlapping between the two scalings: non-local and local (Kolmogorov). This overlapping is based on the very nature 
of the stability exchange between the two statistical regimes. Therefore, it is not a simple task to determine the scale 
r c from the S' n (r)-data. A more fine information one can infer studying 

D LLL (r) = (Au 3 ) 

(cf Eq.(8) and see Q , [H| , [lH ) . Unlike S3 most positive contributions are canceled by negative ones in Dlll, and only 
the slight asymmetry of the Au probability density contributes to Dlll- However, this asymmetry has a fundamental 
nature. Therefore, one can expect that the above used dimensional considerations can be also used for Dlll for the 
both non-local and local regimes. For the local (Kolmogorov) regime one has \l\ . flEI] . fl8l] 

-Dlll ~ E r (25) 

while for the degenerate scaling 

-Dlll- IHr/v)} 5 (26) 

Figure 7 shows —Dlll (as circles) and —Dlll/st (as crosses) calculated using data from the direct numerical 
simulation of homogeneous steady three-dimensional turbulence [1 81 ) . The solid curve in this figure corresponds to the 
best fit by equation (26). In this figure one can see that the overlapping between the two regimes exists even for the 
Dlll (cf Fig. 2). However, a generalized energy input rate related to Dlll could give a key to our problem. 
Let us consider the Navier-Stokes equations for a viscous incompressible fluid, given by 

du t du, I dp d 2 Ul , 

~at + Uj a^ = -pdVi + + /l(x ' t] > (27a) 



where f(x, t) is random force, v is kinematic viscosity and p is the fluid density, f will be assumed to be Gaussian 
with zero mean and a rapidly oscillating character, or a (S-correlation in time. The second-rank correlation tensor 

</i(x + r,t + T)/ J (x,t)> ~ Fij (r)6(T). (28) 

defines such fields. It was shown in [24j using Eqs. (27) and (28) that S2 and Dlll are related by the equation 
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FIG. 7: The DNS data [l8l ] (R\ = 460) are shown against r/rj as the symbols: — D^LL-circles and —Dlll /er-crosses. The 
solid curve is the best fit of (26) corresponding to the degenerated scaling. 




FIG. 8: —Dlll against r/rj. The same data as in Fig. 7 but with higher resolution in a small vicinity of the x m scale. The 
solid curve is the best fit of (26) corresponding to the degenerate scaling. 

dS 2 f r 

Dlll = 6^—^ - — A \ x 4 F u (x)dx, (29) 
dr r 4 J a 

It was also shown in [24| that 



F«(0) = 2(e). (30) 
Fu corresponds to an external energy input rate (cf Novikov's relation [23|: (/i(x, t)vj(x', t)} = | ^ (x — x')). 
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FIG. 9: The same as in Fig. 8 but for R\ = 125. 



For r«£ [24 



D 



LLL 



dS 2 
dr 



When the viscous term can be neglected this gives the Kolmogorov's law (25). 
One can formally rewrite (29) as 



D 



LLL 



= — j / x 4 F(x)dx, 



where the generalized energy input rate is defined as [151 ] 

F{x) = F u {x)- 



3u d(x dS2/dx) 
x 4 dx 



(31) 



(32) 



(33) 



Without the second (viscous) term in the right-hand side of Eq. (33) one obtains the Kolmogorov law (25). 
Therefore, one can associate the first term in the right-hand side of the Eq. (33) with the local (Kolmogorov) regime 
while the second term can be associated with the non-local one. Gradients of these terms can be associated with local 
and non-local interaction strengths respectively. Then balance of the local and non-local interaction strengths should 
be reached at the point given by the equation dF(x)/dx = 0. Let us denote position of the balance point 
One can expect that the crossover scale r c between the local and non-local regimes (see Section IV) coincides with 
the balance point of the local and non-local interaction strengths x m , i.e. r c — x m . It was shown in [151 ] that position 
of the generalized energy input rate maximum x m = r m /1.22, where r m is position of —Dlll(t) maximum. Since 
r m can be obtained from the data this gives us a possibility to estimate r c from the data as well. For instance, we 
know [l5j that r m /rj ~ R^ 7 , hence r c /rj ~ i?°' 7 . In Fig. 7 we indicated position of the scale x m by an arrow and one 
can see that overlapping of the two regimes can be substantial even for Dlll- In figure 8 we show the same data 
as in Fig. 7 but with a considerably higher resolution in a small vicinity of the x m scale. The degenerate scaling 
fit (solid curve) indeed seems to begin its declination from the data in a small vicinity of the x m -sc&\e in this case. 
To support this point we show in figure 9 analogous situation observed for considerably smaller Reynolds number 
Re\ — 125 (see also next Section). 



VII. PASSIVE SCALAR 



The inertial-convective regime of a passive scalar field transport in isotropic turbulence for large Reynolds numbers 
is usually described using Obukhov-Corrsin theory This theory predicts scaling spectrum for the passive scalar 
(6) fluctuations with the scaling exponent equal to "-5/3" 



E 6 {k) 



-1/3 



-5/3 



(34) 



11 




FIG. 10: S® against r/rj. The DNS data (circles) of homogeneous isotropic turbulence described in [28|], Peclet number P\ = 427 
and the Schmidt number is unity (i.e Pa = R\). The solid curve is the best fit to the degenerate scaling (40) (the dashed 
straight line indicates the Obukhov-Corrsin ordinary scaling 



where ee is the average value of dissipation rate of scalar variance, eg (cf Eq. (9)). In the experiments, however, this 
regime is not usually observed (especially, for moderate Reynolds numbers and in turbulent shear flows). In the last 
case the experiments indicate a strong dependence of the passive scalar spectra on the Reynolds number (see, for 
instance, Refs. [25j], [26| and references therein). In the vein of the Section IV the dimensional considerations applied 
to the non-local asymptotic result in the power-law spectrum 

E e (k) ~ — k' 1 (35) 
s 

(cf Eq. (10) and [H| , i H3 ) ■ Accordingly, both scaling regimes (34) and (35) can be corrected (for simplicity we 
restrict ourselves by the case when the Schmidt number is unity [28|]). The first order correction is 

E e (k) ~fc- 5/3 [ln(fc d /fc)] 7 (36) 

and 

£ 9 (fc)~fc- 1 [ln(fc d /fc)f (37) 

for the local and non-local regimes respectively (7 and f3 are dimensionless constants, different from those in Eqs. 
(H),(12)). 

One can estimate the second order structure function for the passive scalar fluctuations 

S 2 V) = (|A0| 2 ) = <|0(r + x)-tf(r)| 2 > (38) 

as 

S 9 2 (r) ~ / E e {k)dk (39) 

Jl/r 

(cf Eq. (13)). Then for r < r c using Eqs. (37) and (39) one obtains 

S$(r) ~ Mr MP (40) 
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FIG. 11: against 7-/77. The DNS data (circles) of homogeneous isotropic turbulence described in [2^], Peclet number P\ = 427 
and the Schmidt number is unity (i.e P\ — R\). The solid curves are the best fit to the degenerate scaling (44). 



with C2 = £ + 1 (cf Section IV). 

Figure 10 shows Sf against r/77 for the DNS data of homogeneous isotropic turbulence described in [28[, Peclet 
number P\ = 427 and the Schmidt number is unity (i.e. P\ = R\). The solid curve is the best fit to the degenerate 
scaling (40) (the dashed straight line indicates the Obukhov-Corrsin ordinary scaling [?j],[28|). As in Section IV we 
can generalize (39) introducing an effective spectrum E„(k) 



>l/r 

(E^ik) = E e {k)). Then using the dimensional considerations we obtain for the non-local regime 



The first order correction is 



Substituting (43) into (41) one obtains 



S«(r) = (|AC> ~ / K(k)dk 

Jl/r 



, ~ v nil 

K(k) ~ (^) 

, P \ n/2 

E e n {k) ~ (-f) fe-X*a/fe)]^ 
#(r) - (f ) Kr/»?)] c - 



(41) 

(42) 
(43) 
(44) 
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FIG. 12: Exponents extracted from Fig. 11 (circles) and for the fully developed atmospheric turbulence (ordinary scaling in 
the inertial interval, triangles .29]).) 

with ( n = Pn + 1. 

Figure 11 shows S„ against r/rj for the DNS data of homogeneous isotropic turbulence described in Peclet 
number Pa = 427 and the Schmidt number is unity (i.e Pa = Pa)- The solid curves are the best fit to the degenerate 
scaling (44). 

One can readily expand the conclusions of the section V regarding Extended Self-Similarity on the case of passive 
scalar and we show effectiveness of this in figure 12. The data for degenerate scaling (circles) were taken from the 
degenerate scaling shown in Fig. 11 whereas the data for the ordinary scaling (triangles) were taken from a high 
Reynolds number atmospheric experiment (29j . 

The overlapping phenomenon (Section VI) takes also place for the 5„ as for S n . In (Ad 3 ) like in Plll most positive 
contributions are canceled by negative ones and only the slight asymmetry of the A9 probability density contributes 
to (AO 3 ). However, similarly to Plll this asymmetry has a fundamental nature. Therefore, one can expect that 
(A9 3 ) also can be used in order to determine r c . Moreover, since (A9 3 ) as function of rjr\ has a turnover point, as one 
can see in figure 13, the turnover point should indicate the point where the actual declination of the degenerate scaling 
fit from the data should begin. In the turnover point derivative of the function (A9 3 ) has its minimum. In the inset to 
Fig. 13 we show with high resolution the derivative in a vicinity of its minimum. We also indicate with arrow the posi- 
tion of local maximum of the generalized energy input rate x m ~ r c calculated using P^ll for this case (cf Section VI). 

One can see that the suggested perturbation theory with the degenerate scaling is an appropriate tool for description 
of data in the near-dissipation range. For modest Reynolds numbers the degenerate scaling applies to nearly the entire 
range of scales. Crossover between the degenerate and ordinary scaling provides an explanation to the Extended 
Self-Similarity. Interplay between local and non-local interactions can be considered as a possible hydrodynamics 
mechanisms of these phenomena. 
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FIG. 13: (A# 3 ) against rjr\. The DNS data (circles) of homogeneous isotropic turbulence described in [28j |. Peclet number 
P\ = 427 and the Schmidt number is unity (i.e Pa = R\)- The solid curve is the best fit to the degenerate scaling. The inset 
shows with high resolution the derivative of (A# 3 ) in a vicinity of its minimum. The arrow indicates position of the generalized 
energy input rate maximum x m ^ r c calculated using Dlll- 
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